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Abstract 

New CP violating physics in b —► s transitions will modify the CP asymmetries in B decays into 
final CP eigenstates {fKg, rj'Kg, ^°Kg, aiKg, p°Kg and yKs) from their Standard Model values. 
In a model independent analysis, the pattern of deviations can be used to probe which Wilson 
coefficients get a significant contribution from the new physics. We demonstrate this idea using 
several well-motivated models of new physics, and apply it to current data. 
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I. INTRODUCTION 


CP asymmetries in the decays of neutral B mesons into final CP eigenstates / exhibit a 
time-dependent behavior (for a review of CP violation in meson decays, see 1]), 

Af(t) = Sfsm(Amt) — CfCos(Amt). (1) 

The Standard Model predicts that for most of the decays that proceed via the quark tran¬ 
sitions b —> q'q's ( q' = c,s,d,u ), the following relations hold to a good approximation: 


—r]fSf~sm2/3, Cf ~ 0, (2) 

where r]f = ±1 is the CP eigenvalue for the final state /, and /3 = (f>\ = 
arg [— {V c dV* b )/(VtdVtb)\- New physics effects can appear in two ways. First, new physics 
in 'H/\b =2 contributes to the B° — ~B° mixing amplitude. Such a contribution shifts all S/ s 
in a universal way: The Sf asymmetries remain equal to each other, though different from 
sin 2/3. The C/’s still vanish. Second, new physics in Hab=i can contribute to the b —> qqs 
transitions (q = s, d, u ) Q (The tree leve, tuition 5 -e M * unlMy to get a eigaiflcaat 
contribution from new physics, and consequently the asymmetries S^k and C^k would not 
be modihed by 'Hab=i-) Such a contribution would lead to interesting consequences: 

1. The S'/’s would be different from each other and from S^k- 

2. The C/’s would be different from each other and from zero. 

Thus, in the presence of new physics, a pattern of deviations, SSf = —r]fSf — S^k 7 ^ 0 and 
Cf 7 ^ 0, will arise. 

Given a set of experimental ranges for <5 Sf and C/, one would like to interpret these data 
in terms of new physics. There are two ways to proceed. First, one can choose a model of 
new physics and analyze whether the model can accommodate the data. The second way 
is model independent. The effects of new physics can be described by a modification of the 
Wilson coefficients in the operator product expansion for A B = 1 interactions. Thus, one 
can fit a set of Wilson coefficients to the data and learn which operators can account for the 
observed deviations. In this paper, we study this second, model independent, method. 

While our analysis has the advantage of being model independent, it suffers from two 
limitations: 
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TABLE I: World averages of the experimental results for the CP asymmetries in neutral B decays 
via the quark transitions b —* q'q's modes, q' = c, s, d, u. 


/ 

-r/cpSf 

Cf 

Refs. 

i’Ks 

+0.73 ± 0.04 

+0.03 + 0.03 

[6 ,H 

<PK S 

+0.35 + 0.21 

-0.04 + 0.17 

[8,9] 

rfK s 

+0.43 + 0.11 

-0.04 ± 0.08 

[a, to] 

ir°Ks 

+0.34 + 0.29 

+0.08 + 0.14 


uK s 

+0.56 + 0.32 

-0.49 + 0.25 

[9,12] 


• The number of Wilson coefficients is larger than the number of measured CP asym¬ 
metries. Consequently, an analysis that is entirely generic is impossible to carry out 
at present. Hence, one has to assume that only a subset of all possible operators 
are modified. In this work we consider some simple cases, where the new physics is 
parameterized by a single complex parameter. The cases that we study are motivated 
by actual models of new physics. Extensions to a larger number of new physics pa¬ 
rameters is left for future work. It is also possible to include in the analysis additional 
data, beyond the CP asymmetries in decays into final CP eigenstates. This extension 
is also left for future work. 


• To interpret the data in terms of Wilson coefficients, one has to know the values of the 
(mode-dependent) hadronic matrix elements of the operators. At present, there is no 
first principle calculation of the matrix elements that has been tested to a high level 
of precision. Thus, hadronic uncertainties prevent a clean theoretical interpretation. 
We use factorization ^,0,1^] for our analysis. 1 

A number of relevant measurements exists already. The experimental situation is sum¬ 
marized in Table |U 


1 In this work we perform our calculations to leading order in a s and drop 1 /mb corrections, except for the 
chirally enhanced terms related to ag and cig. In this approximation, QCD- and naive-factorization are 
identical. We note, however, differences in the expressions for the / = r/'^K modes between refs. dJ and 
fs]: A term oc /“ {() (<26 — ^ag) should be omitted from the first, while a term cx /))(,) (<27 — ag ) should be 
added to the second. 
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The plan of this paper is as follows. In section UTl we introduce our formalism and evaluate, 
to leading log approximation and within the QCD factorization approach, the Standard 
Model (SM) predictions for the CP asymmetries. In section Hill we analyze the effects of new 
physics. We start with a generic, model independent analysis, using the operator product 
expansion for A B = 1 operators. We then focus on scenarios where the effects of new physics 
depend on a single complex parameter. We give three explicit examples of such scenarios. 
In section HVl we apply our approach to present data. We conclude in section lYl 


II. FORMALISM AND STANDARD MODEL PREDICTIONS 


f Q 


We follow the notations of ref. |3[. We consider the following effective Hamiltonian for 


A B = 

= ±1 decays: 






7 " 0 fr — 

g f 

V2 

E 'y.us (cio?+C205 • 

10 \ 

T 'y ' CiOi + (AyO7~ f + CsgOsg | + h.c., 

i= 3 / 

( 3 ) 

with 










Oi = ( pb) v _ A (sp)v _ A , 

02 = 

(. P/ 3 b a )v-A(s a Pi 3 )v-A, 




o 3 

= (sb) v _ A ^2(qq)v_ A , 

0 4 = 

{s a bp)v-A y 1 (g/?gaV-A, 




o 5 = 

= {sb)v-A^(qq)v+A, 

q 

o 6 = 

(s a bp)v-A y ^i.qpqa)v+A, 

q 



O7 

_ 3 
2 

(■ sb) V -A^e q (qq) v+A , 

q 

o 8 = 

3 _ 

-{s a b/ 3 ) v _ A '^2e q (q 0 q a )v+A, 

q 



o 9 

3 

= 2 

{sb)v_A^e q (qq) v _ A , 

q 

0 10 = 

q 



O77 

= 

■^^(1 + 75)^, 

0 8g = 

9s m b _ 1 L 

g 2 sa (■*■ + 7-5 jG^b, 

( 4 ) 

where 

; ( QlQ2)v±A : 

= gi 7 M (l ± 75)92, the sum 

is over active quarks, with e q denoting 

their 


electric charge in fractions of \e\ and a,/3 are color indices. 

The CP asymmetries in B —> f decays are calculated as follows. One defines a complex 
quantity A/, 

A , = e-**‘(3 f /A,), (5) 

where 0# is the phase of M 12 , the B° — B° mixing amplitude, and Af ( Af ) is the decay 
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TABLE II: The anomalous dimension matrix at leading log. The rows and columns correspond to 
the twelve operators Oi,...,io, 77 , 8 g defined in eq. 0). 


-2 

6 

2 

9 

2 

3 

2 

9 

2 

3 

0 

0 

0 

0 

416 

81 

70 

27 

6 

-2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

22 

9 

22 

3 

4 

9 

4 

3 

0 

0 

0 

0 

464 

81 

545 

27 

0 

0 

44 

9 

4 

3 

10 

9 

10 

3 

0 

0 

0 

0 

136 

81 

512 

27 

0 

0 

0 

0 

2 

-6 

0 

0 

0 

0 

32 

9 

59 

3 

0 

0 

10 

9 

10 

3 

10 

9 

38 

3 

0 

0 

0 

0 

296 

81 

703 

27 

0 

0 

0 

0 

0 

0 

2 

-6 

0 

0 

16 

9 

5 

6 

0 

0 

1 

9 

1 

3 

1 

9 

1 

3 

0 

-16 

0 

0 

1196 

81 

11 

54 

0 

0 

2 

9 

2 

3 

2 

9 

2 

3 

0 

0 

-2 

6 

232 

81 

59 

54 

0 

0 

1 

9 

1 

3 

1 

9 

1 

3 

0 

0 

6 

-2 

1180 

81 

46 

27 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

32 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

32 

9 

28 

3 


amplitude for B°(B ) —> /. We have 

2Im(A / ) 1-IAjl 2 

’ 1 + |A/| 2 ’ ' 1 + IA,! 2 ' 

The decay amplitudes can be calculated from the effective Hamiltonian of eq. 0 

A, = A, = (/|H«*|B°). 


( 6 ) 

(P 


The electroweak model determines the Wilson coefficients while QCD (or, more practi¬ 
cally, a calculational method such as QCD factorization) determines the matrix elements 

(f\Oi\B°(B°)). 

We perform our calculations to leading-log approximation. In particular, to run the 



Within the SM, we have the following set of Wilson coefficients at leading order: 


Cf M (m w ) = 1, Cf™(m w ) = 0. (8) 
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(Strictly speaking, C^{mw) and C^{rriw) are also different from zero. However, their 
contributions to the decay processes of interest occur at next-to-leading order which we 
neglect in the present work.) The SM contribution to the decay amplitudes, related to 
b —> q'q's transitions, can always be written as a sum of two terms, + AJ, with 

Aj oc V* b V cs and Af oc V* b V us . Defining the ratio aj = e~ n (Af / A^) , where 7 = 0 3 = 
arg [-(V ud V* b )/(V cd V; b )\, we have 

Af 1 = A c f (1 + a u f e^) . ( 9 ) 


For B —* ipKs , the aj term can be safely neglected (its effects are below the percent level) 
and, consequently, 

( 10 ) 

For charmless modes, the effects of the a 1 } terms (often called ‘the SM pollution’) are at 
least of order \(V ub V* s )/(V cb V c *)\ a few percent. They can lead to a deviation of Sf from 
S^pK and of Cf from zero. In Table Hill wejjfiye the values of the a'j parameter (obtained, as 


explained above, by using factorization 


000 ) 


for all relevant modes. Since we perform our 


calculations to leading log approximation, we neglect, for consistency, also power corrections 
within the QCD factorization approach, except for the dominant chirally-enhanced ones 
related to Qq and ag. This means that the strong phases vanish, and that Cf — 0 for all / in 


this approximation. The SM predictions for Sf (with S^k s = 0-73 and 7 = 62° 0 taken 
as input) are also given in Table ITTT1 The values of the input parameters used are specified 
in Table m To estimate our uncertainties, we give the values for n — m& = 4.2 GeV and 
for a varied renormalization scale of p = m b / 2 and /i = 2ro;,. Note that the coefficients 
of the chirally enchanced terms, defined as r x = 2nr? K j and = h s s /(fgm b m s ) for 
the singlets S = 7 , 7 / (see J] for details), evolve with the renormalization scale through the 
running (MS-bar) quark masses. 

An examination of Table II I II shows that the SM pollution is small (that is, at the naively 
expected level of \(V u bV* s ) / {V c bV* s )\ ~ a few percent) for / = (j)Ks, r/Ks and 7r °Ks- It 
is larger for / = rjKs, ujKs an d P°Ks- I 11 these modes, a^ is enhanced because, within 
the QCD factorization approach, there is an accidental cancellation between the leading 
contributions to Aj. The reason for the suppression of the leading Aj piece in / = pK , uK 
versus / = 7 t°A' is that the dominant QCD-penguin coefficients a 4 and a 6 appear in A c ^ p ^ K 
as (a 4 — r x a 6 ) and in A° 0R as (a 4 + r x a 6 ). Since r x ~ 1 and, within the Standard Model, 
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TABLE III: The parameters, calculated in QCD factorization at leading log and to zeroth order 
in A /m b (except for chirally enhanced corrections), and the SM values of Sf for p = m b and in 
parentheses the respective values for p = 2 m b (first) and p = m b /2 (second) if different from the 
central one. 


/ 

a u f 

-rjcpSf 

i/K s 

0 

0.73 

ct>K s 

0.019 

0.75 

tt°Ks 

0.052 [0.094,0.021] 

0.79 [0.83,0.76] 

t]K s 

0.08 [0.16,0.02] 

0.82 [0.88,0.76] 

r/Ks 

0.007 [-0.006,0.019] 

0.74 [0.72,0.75] 

uK s 

0.22 [0.37,0.04] 

0.92 [0.98,0.78] 

p°K s 

-0.16 [-0.32,0.005] 

0.49 [0.19,0.74] 


TABLE IV: Input parameters used in the calculations. 


Form Factors 

F B ^ n ( 0) 

0.28 

F b ^ k { 0) 

0.34 

F b ^( 0) 

0.23 

F b ^\0 ) 

0.19 

A b ^( 0) 

0.33 

A b ^p{ 0) 

0.37 


Decay constants (MeV) 

fn 

131 

Ik 

160 

Jv 

108 

fS 

-111 

$ 

89 

Fr,' 

136 

F 

187 

fp 

209 


r x -factors, quark masses 

r x (m b ) 

1.170 

r X { m b) 

1.230 

r x( m b) 

1.241 

m b (m b ) 

4.2 GeV 

m s { 2GeV) 

110 MeV 

(m u + m (i )(2GeV) 

9.1 MeV 


CZ 4 ~ a 6 , there is a cancellation in A C ,^ K while there isn’t one in A c / 0k . The suppression for 
A^ k with respect to A C , K has a different reason: it is due to the octet-singlet mixing, which 
causes destructive (constructive) interference in the r)(r)')K penguin amplitude jl3|. 

We stress that the numerical results presented above are often sensitive to the approxi¬ 
mations that we make and to the values of the input parameters. We will refine them in the 
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future by going beyond the leading log approximation and taking into account uncertainties 
in the input parameters other than the scale fi. At present they should be taken as indi¬ 
cation to how close to S$k s one should expect the various S/s to be, but not as accurate 
predictions. Our findings are compatible with the next to leading order results in the SM 


given by 


Q 


We conclude this section by adding some general considerations concerning the accuracy 
of our approximation and the stability of our results. We discuss two main points. First, 
the validity of our approximation for branching ratios and direct CP asymmetries. Second, 
the presence of large power corrections in our analysis. 

Regarding the first issue, there is a substantial effect from NLO corrections in factorization 
as far as the branching ratios are concerned. However, the impact of such NLO corrections 
is very moderate for the asymmetries Sf, which are the main focus of our study. This 
can be demonstrated explicitly in the SM b y co mparing our results (Table IIIII) with the 
more complete, full-fledged NLO analysis of |Ij|. In fact, even the estimate of theoretical 
uncertainties, which we obtained by a variation of the renormalization scale, is in good 


agreement with Q, w.rere also other sonrees ot nneertamty are me.uded. This clearly 
demonstrates that our approach gives a reasonable picture for the asymmetries Sf, including 
the issue of uncertainties and the stability against subleading effects. 

In addition, our framework is still consistent with the observed direct CP asymmetries 
Cf, which are compatible with zero. This is the case in general for direct CP asymmetries in 
all B decays, where Acp(K + n~) is so far the only exception. Even in this case the observed 
direct CP asymmetry is not very large (about 10%). 

Regarding the second issue, we stress again in this context that the dominant power 
corrections related to Qq and Qs are included in our analysis (as already mentioned in 
footnote 1). Further power corrections, in particular those from annihilation topologies 




can be estimated U to be typically of the order 10 - 20% with respect to the faetorizable 
prediction. Hence the impact of such effects is subdominant. This is also confirmed by the 
comparison with mentioned above. 



III. NEW PHYSICS 


A. Generic analysis 


For the purpose of discussing new physics, it is convenient to define a phase /3 e g: 


A 


— _ e -2if3eff 


ipK s = 


( 11 ) 


In writing down eq. m we are making the very plausible assumption that the quark 
transition b —> ccs is dominated by a single weak phase. This is clearly the case within 
the SM, where the phase of V cb V* s dominates to better than one percent. But even with 
new physics we do not expect significant new contributions to SM tree level processes. The 
approximation of a single weak phase may be not as good as in the SM, but it is still expected 
to be a very good approximation. On the other hand, eq. m allows for new physics in 
B° — B° mixing, in which case (3 e g- ^ j3. 

Contributions of new physics to H c s of eq. © modify the Wilson coefficients C t . (New 
physics can also introduce additional operators. We do not consider this possibility here, but 
our analysis can be generalized to this case in a straightforward manner.) This modification 
can be parametrized as follows: 

Ci(m w ) = Cf M (m w ) + Sie l6 \ (12) 


where Cf M is the Standard Model value of the Wilson coefficient, is real and positive, and 
9i is a phase in the range [0, 27 t], 

The modification of Aj from the SM expression [eq. ©] due to the new physics contri¬ 
bution [eq. (HI] can always be written as follows: 


A f = A 


f 


1 + a u f e + 


E 


+ ifce*) e<e 


i6i 


(13) 


n Ta 


Die 0 we give the values of the parameters within the QCD factorization approach 
5|], The Iff t for % — 3,..., 10 are related to the b c fi through b v j t = {\V ub V* s \ / \V cb V* s \) . 

Note that further b < j 1 = 1 and b v j x = aj. 

With new physics, the deviation of Sf of charmless modes from S^k depends on b r j - l , 6^, 
e t and 9i. If enough measurements of Sf and Cf asymmetries become available, one will 
be able to find the values of the £. t - and 6 t - parameters that account for the data. Given a 


9 







small set of measurements, one can still find £* and dj if only a small number of operators is 
affected by the new physics in a significant way. In the next section, we focus on a simple 
(but well motivated) case in which the contribution of the new physics can be parametrized 
in terms of a single complex parameter. 

B. A single complex parameter 

Suppose that all the parameters of eq. m are related to a single complex param¬ 
eter, that is, 

£ie ldi = Xi£e %6 for all i, (14) 

where, in a given model of new physics, the Xj-constants are computable (see section IIII Cl 
for three specific examples). In such cases, eq. (11311 takes the following form: 

A f = A) [ 1 + ay' 1 + (b c f + by" 1 ) £e ie ] . (15) 

The analysis is simpler if \b c j\ \bj\. This is the case for all generic combinations of operators 
Oi in eq. (J3|) that do not involve i = 1,2. We proceed therefore in discussing the case where 
a combination of operators within the set i — 3,..., 10, 7y, 8 g is affected by new physics. 
For two classes of interesting scenarios, simple analytical expressions can be obtained for 
the modification of Sf and Cf. 

(i) The new physics contribution is dominant. Explicitly, 

e\b C f \ 3> \a U f\, 1 . (16) 

The shift in all modes where eq. m holds is universal and depends only on 6\ 

— r]fSf ~ sin 2/? e ff cos 29 + cos 2/3 e s sin 2 d, 

c, ^ o. (17) 

(ii) The new physics contribution and the SM pollution are small. Explicitly, 

e\b C f \ -C 1 and \a U f\ -C 1 . (18) 

The shift is mode dependent and depends on esind: 

— rjfSf ~ sin2/3 eff + 2 cos2/3 eff [' TZe(b < j)£sm9 + 7 Ze(a^) sin 7 ] , 

Cf ~ — 2Im(b C f)£ sin d — 2Im(a U f) siny. (19) 
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Since the b c j values are known (see Table 0, the pattern of deviations is predicted and can 
be used, by comparing with the data, to close in on the set of operators that is responsible 
for the deviations and to extract £ sin 9. 

In case that none of these approximations applies, one can still find the pattern of devia¬ 
tions by numerical evaluation for any values for £ and 6. By comparing these theoretical SSf 
and Cf to the data, allowed regions in the e — 9 plane (for each scenario) can be determined, 
and a best fit point (among all scenarios) can be found. We demonstrate this procedure in 
the next section by using present data. 


C. Specific Scenarios 

We now consider three explicit examples of new physics models. All three examples are 
well motivated and fulfill eq. da- „ 

(i) The dominant effect of the new physics is through Z-penguins [12]. Then, the Wilson 
coefficients are modified as follows: 


C^imw) = ^ z e id % 

o 


C^ F (m w ) = — sin 2 0 W £ z e l9z , 
Cr(m w ) = -|(1 - sin 2 9 w )£ z e l6z 


Here 


e,e‘ e ‘ = 3 


J sb 


( 20 ) 


( 21 ) 


4^v a v,;' 

where Z s b parametrizes the coupling of the Z° boson to the left handed strange and bottom 
quarks. The consistency of the measured BR(6 —> s£ + £~) with the SM prediction requires 
that \Z s b\ < 0.08 j3] which gives 

£,< 0 . 02 . ( 22 ) 


(ii) Kaluza-Klein excitations of gluons in RS models [2P] couple mainly to the third 
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generation quarks and contribute to all gluonic penguin operators: 2 

C± P (m w ) = C^ p (m w ) = £ k e ldk , 

C^ p (m w ) = Cf p {m w ) = ~£ k e ldk . (23) 

Here 

E ‘ e "‘ = -™ JMg) (sb) iffe* (24) 

where D sb is the rotation matrix from interaction to mass basis for the left-handed down 
quarks, x = ®(X) is a model dependent parameter (related to bulk fermion masses), and 

Mq is the mass of the lightest excitation. With M G > 1 TeV and \D sb \ ~ \ V t i,V t *\, we expect 


0.02. (25) 

(iii) Enhanced chromomagnetic operator 0 can be parametrized as follows: 


Csf(m w ) = £ g e v 


(26) 


(We still make the approximation that Cf^(mw) = 0 and neglect further Cf^m?,) induced 
by 0\.) From BR(R —■> X s g) < 0.09 22j, 23], we have 


£ g ~ 


(27) 


Note that in this case the new physics contribution enters at next to leading order. Yet, the 
modifications of the CP asymmetries can be substantial. An experimental improvement by a 
factor of a few in the upper bound on BR(H —► X s g) will restrict the potential modifications 
in this scenario in a significant way. 

We have calculated the b c f and 67 parameters of eq. (fT51) for each of the three scenarios. 

rim 

We use the factorized decay amplitudes from [3|, [4|, |5|. (We use asymptotic distribution 
amplitudes in the analysis of the scenario with an enhanced chromomagnetic operator.) The 
results are presented in Table IVI1 The Uj parameters are quite stable under variation of the 
scale. As anticipated in section lTll Bl in all three scenarios we have bj = (\V ub V* s \/\V cb V c * |) bj-. 


2 In ref. fell ] it is argued that, quite generically in the RS framework, the set of operators in eq. m is 
subdominant and constrained to give only a small contribution to the decays in question. However, this 
issue is model dependent, and eq. also represents a viable scenario. 
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Given the upper bounds on e for the three scenarios - £ z < 0.02, e k < 0.02 and e g < 1 
- we learn that in none of the three examples can the new physics strongly dominate over 
the SM contributions, that is, eb c j a'j, 1. Thus, if the data makes a convincing case for a 
universal pattern of deviations, it would mean that either this pattern is accidental or that 
the new physics is different from the three specific scenarios discussed here. 

On the other hand, the interpretation of the results presented in Table I VI I is straight¬ 
forward in the other limit discussed in the previous subsection. For the four modes 
/ = 7T °K S , VK s , ri'K s and 4>K S , we have, say, \b c f e\ < 0.3 for £ 2 < 0.01, e k < 0.005, 
and e g < 0.17. In any of these cases, the deviation from the SM prediction is proportional 
to b c r [see eq. dH ]. Consequently, there is a distinct pattern of deviations for each of the 
three scenarios. For example, the deviation of S v >k s from the SM prediction is a factor of 
four smaller than that of S v k s in the first scenario (£ z ^ 0), a factor ~ 1.3 smaller in the 
second (£ k ^ 0) and a factor of ~ 1.3 larger in the third (£ g ^ 0). The deviation of S^Ks is 
similar or somewhat smaller than that of S n K s hi the first two scenarios, but a factor of ~ 1.9 
larger in the third. The second scenario has the characteristic feature that the SM deviations 
of Srf,K s and S^k s are very similar, whereas in scenario i (iii) S^Ks is a factor of 3 (1.6) larger 
than Sg’Ks- It is interesting to note that in all three new physics examples, the deviations 
for b c j£ -C 1 in the modes that have been most accurately measured, / = tt°Ks, ri'Ks and 
4>Ks, and also r/Ks, are in the same direction, that is, the four Sf asymmetries are either 
all larger or all smaller [depending on sign(sin6 ) )] than the SM prediction. 

To demonstrate how the pattern of deviations probes new physics, we perform the fol¬ 
lowing exercise. For each of the three scenarios discussed in this section, we take £; at the 
phenomenological upper bound, that is = 0.02, £ k = 0.02 and £ g = 1. Since the new 
physics contribution is now as large as it can get, the patterns of deviations are expected to 
be more significant. We then pick the two values of 6i for which S n R s — 0.34, the experi¬ 
mental central value. (There is no particular reason for selecting S n K s - O ur aim is just to 
demonstrate the idea that the pattern of deviations is sensitive to new physics.) We present 
the results in Fig. |H One can clearly see that the six patterns are different from each other, 
so that each scenario can in principle be easily tested. 
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FIG. 1: Pattern of Sf asymmetries for the three scenarios presented in section rrrm with fixed 
SttK s = 0.34. The symbols represent • = ir, ■ = r), A = ?/, T = <f>, X = u and + = p. The upmost 
pattern gives the asymmetries in the SM for ^ = mb/2,nib and 2mb- The other six patterns, from 
top to bottom, correspond to ( £ z ,d z ) = (0.02,0.85), (e 2 ,6b) = (0.02,1.79), (e^,^) = (0.02,0.12), 
(sfci^fc) = (0.02,2.52), ( £g,0 g ) = (1,3.27), and ( £ g ,0 g ) = (1,5.64). The vertical bar represents the 
experimental value for S^k s ■ 

IV. PRESENT DATA 

At present, the B factories provide constraints on eight relevant CP asymmetries. The 
data are presented in Table [U We would like to test various scenarios of new physics by 
performing a \ 2 fit to the data. Since we are mainly interested in demonstrating the potential 
power of the data to probe new physics in the future, we make several simplifications: 

1. We work within the framework defined in section HUB, that is, new physics contribu¬ 
tions that depend on only a single complex parameter. 

2. Since the measurement of S^k s is very accurate, we fix the value of sin 2 f3 e s to its 
central value, and do not include it as a fit parameter. Furthermore, since experimental 
data Q disfavor cos 2/3 eff < 0, we allow only positive cos2/3 eff values. 
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We thus calculate, for various scenarios, 


x 2 M) 


E 


Xf p - xf F (e, 9) 


a 


exp 


(28) 


where X c ^ p and (7j Xp are, respectively, the central value and the one sigma range of the 
experimental measurements of Sf and Cf for / = tt°Ks, r/'Ks , 4>Ks and cuKs, and Y^ F (£, 9) 
are the (central) theoretical values of Sf and Cf calculated within the QCD factorization 
approach for given values of e and 6. Only experimental uncertainties are taken into account. 
Since there are eight observables and two free parameters, y 2 has a six degrees of freedom 
distribution. 

We consider eleven cases. Firstly, we investigate the possibility that the new physics 
contributes significantly via only a single Oj, where i = 3,..., 10. (Specific models of new 
physics do not provide a strong motivation for such scenarios. Yet, this investigation is 
useful in the two limits discussed in section HIT Bl If the new physics is dominant and 
depends on a single weak phase, we can see the universality of its effects [eq. (IT71) ]. If the 
new physics contribution is small, the effects are well approximated by a sum over the shifts 
from the various operators [eq. (|19|) ].) In each case, we perform our analysis for e < 5 and all 
9 e [0, 27 t]. Secondly, we investigate the three specific examples presented in section UlI Cl 
We do this for £ within the phenomenologically allowed ranges and for all 9 values. 

We use the y 2 fit in two ways. For each case we find the favored regions in the e — 9 
plane. The results are shown in Figs. El and 0 Since the fits to individual operators O 3...10 
give very similar results for large e, we only show the ones for O 4 , O 5 and Og. Then we test 
the overall goodness of the fit for each case and learn whether the possibility that the data 
are accounted for by a new physics contribution to this (set of) operator(s) is disfavored. 
We draw the following conclusions: 


1. For almost all the Wilson coefficients, a very good fit is obtained for a high e value, 
such that eq. (11 (ill holds and the shift is universal [eq. ()17jl ]. The reason is that the 
central values of the three best measured Sf asymmetries (/ = rj'Kg, 4>Ks, vr °Ks) show 
an approximately universal pattern: 5Sf ~ —0.4. 

2. For operators where eq. m gives the best fit, the preferred value of 9 corre¬ 
sponds to sin 2(/3 + 9) ~ 0.4. Note that a four-fold ambiguity arises, giving 9i ~ 
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55°, 170°, 235°, 350°, which can be seen in Fig. |2] Had we allowed for cos2/3 e ff < 0, 
there would be eight possible solutions for high e. 

3. For each of the three specific scenarios that we considered, there exists a region in the 
e — 6 plane where a good fit is obtained. Note that within the operator basis we use, 
in the first two scenarios with modified four-Fermi operators a phase smaller than n 
is preferred, whereas n < 6 g < 2n is favored. 

4. Naively, the probability that there is no new physics contribution affecting the mea¬ 
sured S/s is small. (The SM point, £ = 0, lies in the white area.) One has to bear 
in mind, however, that we used a rather crude approximation, and that there are 
large uncertainties in our evaluation of the matrix elements. Furthermore, the dis¬ 
crepancy between S^ s = 0.43 ± 0.11 and S^ <s = 0.74 plays a significant role in 
this result. Note however that the Belle and Babar results for this observable are not 
quite consistent. Inflating the error according to the PDG prescription would yield 
S^/Ks = 0-43 ± 0-17 and a better SM fit. (Now the SM point lies in the light grey 
area.) 
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FIG. 2 : Allowed region in the e, — 9 t plane for new physics in a single operator: (a) O 4 , (b) O 5 and 
(c) O 9 . The black, dark grey, and light grey regions correspond to probability higher than 0.32, 
0.046, and 0.0027, respectively. 
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TABLE V: The b c ^ {i = 1,..., 10) and 2 parameters calculated in QCD factorization to leading 
log approximation at /i = mj and to zeroth order in A/mj (except for chirally enhanced corrections). 


/ 

7T °K 

r]K 

rf K 

<t>K 

loK 

p°I< 

b fl 

1 

1 

1 

1 

1 

1 

h u 

1 

0.052 

0.08 

0.0072 

0.019 

0.22 

-0.16 

h c 
°f 2 

-0.12 

-0.11 

-0.13 

-0.13 

-0.096 

-0.093 

b h 

0.32 

0.59 

-0.12 

-0.0024 

1.9 

-1.8 

b h 

-0.34 

17 

-26 

-46 

210 

9.3 

b % 

-12 

-7.4 

-17 

-43 

110 

81 

b h 

-12 

-28 

12 

-43 

140 

-55 

b % 

-31 

-39 

-21 

-20 

-51 

-160 

b h 

-19 

-32 

3.8 

22 

80 

-110 

h c 

8 

3.6 

-5.1 

20 

14 

120 

-1.6 

h c 

9 

25 

37 

3.5 

23 

45 

-140 

h c 

°f io 

12 

11 

11 

24 

-36 

-51 






-3 0 . 
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FIG. 3: Allowed regions in the e — 6 plane for the three specific scenarios, (a) enhanced Z-penguins, 
(b) KK gluons in RS models, (c) enhanced chromomagnetic operator. The black, dark grey, and 
light grey regions correspond to probability higher than 0.32, 0.046, and 0.0027, respectively. 
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V. CONCLUSIONS 


5-meson decays that proceed via the quark transition b —> qqs , for q = s,d, u , open a 
window to new physics because the SM contribution is either loop- or CKM-suppressed. In 
particular, the CP violating asymmetries Sf and Cf in these modes may reveal the existence 
of new sources of CP violation beyond the KM phase. 

With new physics, these decays may get significant contributions that depend on a phase 
that is different from the SM phase. Then, not only will these asymmetries be different 
from the SM prediction but also, in general, they will differ from each other. The pattern of 
deviations, SSf — —rjfSf — S^k and Cf, allows us to probe in detail the nature of the new 
physics that accounts for the effect. 

Factorization schemes predict that there is an accidental cancellation between the leading 
Standard Model contributions to the 77 Kg, ojKg and pKg modes (see section HU). The 
suppression of the Standard Model amplitudes makes the corresponding CP asymmetries 
more sensitive to New Physics. At the same time, the sensitivity to hadronic uncertainties 
is also stronger. 

We presented a way to analyze the data in a (new physics) model independent way. In 
a low energy effective theory, the effects of the new physics are expressed as modifications 
of Wilson coefficients in the A B = 1 effective Hamiltonian. The data can be used to find 
which Wilson coefficients are modified and by what (complex) factor. 

The analysis is simplified when the modification of the Wilson coefficients depends on a 
single complex parameter, as is often the case in specific models. We gave analytic expres¬ 
sions for the modifications in some interesting limits. In particular, for C^. ^io,8g> when the 
new physics dominates, a universal shift of the asymmetries arises [see eq. (ED]. If the new 
physics contribution is small, then the deviations (for / = rj'Kg, 4>Kg, tt°K s , TjKg) have a 
pattern that depends in a simple way on the ^-parameters [see eq. pp]. 

We applied this method to present data, using lowest order QCD factorization to calculate 
the hadronic matrix elements. We find that, since the best three measured S'/’s have similar 
central values (~ 0.4) which are, however, different from S^k, an d all measured Cf 's are 
consistent with zero, a good fit can be achieved for almost all operators if they are enhanced 
by the new physics to a level where they dominate over the SM contribution. In three 
specific scenarios that we considered, such new physics dominance cannot be realized in 
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Nature because of phenomenological constraints. Still, a good fit can be obtained for each 
of these scenarios. 

Our analysis can be improved in several ways. In particular, the calculation should go 
beyond the leading log approximation and incorporate power corrections. The combination 
of improved calculations, more accurate measurements and additional modes explored ex¬ 
perimentally, is likely to make the pattern of deviations from the SM (or their absence) a 
powerful probe of flavor and CP violating new physics. 
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TABLE VI: The b c j parameters calculated in QCD factorization to leading log approximation and 
to zeroth order in A/m& (chirally enhanced corrections are included) in the three scenarios: (i) 
Z-penguins [eq. (120J) ]. denoted with subindex z\ (ii) KK-gluons [eq. (123J) ]. denoted with subindex 
k\ (iii) Chromomagnetic operator [eq. ( 1261 ) ]. denoted with subindex g. The first row in each case 
is obtained for g = mb, whereas the second row gives the values for g = mb/ 2 (first) and 2 mb 
(second). 


/ 

7 T°I< 

7)K 

g'K 

<j)K 

ljI< 

p°K 

h c 

-16 

-21 

-5.4 

-16 

24 

57 


[-12,-22] 

[-16,-30] 

[-4.0,-7.6] 

[-10,-24] 

[29, 28] 

[60,65] 

h c 
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-34 

-33 

-61 
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[-33,-59] 

[-22,-51] 

[-24, -48] 
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h c 

°fg 
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1.4 
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